Nonlocal effects on magnetism in the diluted magnetic semiconductor Gai-j^Mn^^As 
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The magnetic properties of the diluted magnetic semiconductor Gai_i,Mna;As are studied within 
the dynamical cluster approximation. We use the k ■ p Hamiltonian to describe the electronic struc- 
ture of GaAs with spin-orbit coupling and strain effects. We show that nonlocal effects are essential 
for explaining the experimentally observed transition temperature and saturation magnetization. 
We also demonstrate that the cluster anisotropy is very strong and induces rotational frustration 
and a cube-edge direction magnetic anisotropy at low temperature. With this, we explain the 
temperature-driven spin reorientation in this system. 

PACS numbers: 75.50.Pp, 75.30.Gw, 78.55.Cr 



The discovery of high temperature ferromagnetism in 
diluted magnetic semiconductors (DMS) has stimulated a 
great deal of attention [ij . The interest in these materials 
is due to possible applications in spintronics as the 
source of a spin polarized current or as the base material 
for a chip that can simultaneously store and process data. 

In spite of extensive studies, our understanding of fer- 
romagnetism in these systems is far from complete Q. 
There are a few serious difficulties in the theoretical 
study of DMS: (i) The magnetic interaction between lo- 
cal magnetic moments and itinerant carrier spin, which 
is responsible for the high transition temperature (Tj,), is 
strong and outside the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) regime, (ii) There exists strong disorder from 
the random distribution of magnetic ions, (iii) Nonlo- 
cal effects are expected to be crucial judging from the 
spatially oscillating and anisotropic magnetic interaction 
predicted in theory [3-01 and observed in experiments Q . 

The mean-field study by Dietl et al. [9^ captures the 
main features of DMS systems qualitatively and some 
even quantitatively. However, it ignores strong corre- 
lations, disorder effects, and spatial fluctuations, and 
fails to describe some DMS materials, such that subse- 
quent studies have brought their approach into question 
[lol [in . Studies 01 based on the dynamical mean-field 
theory (DMFT) [13l| have made considerable improve- 
ments by including strong correlation and disorder ef- 
fects. However, the local nature of the DMFT presents 
severe limitations when studying this system. The effects 
of short-range fiuctuations and spatial correlations were 
shown to be important in the classical Heisenberg model 
[14]. In this letter, we show that nonlocal effects may be 
equally important for the itinerant carriers, which me- 
diate the effective interaction between local moments in 
DMS. 

The dynamical cluster approximation (DCA) [ll] sys- 
tematically incorporates nonlocal effects as the cluster 
size {Nc) increases while retaining strong correlations. 



When Nc = 1, the DCA is equivalent to the DMFT, 
and exact results are approached as Nc ^ oo. Since all 
the possible disorder configurations are considered in a 
cluster, the DCA is also a better approximation for disor- 
der average than the coherent potential approximation or 
DMFT by including multi- impurity scattering terms [l6| . 
Thus, the DCA is an ideal method for studying DMS sys- 
tems. In this Letter, we study the magnetic properties 
of the prototypical DMS system Gai-xMn^jAs using the 
DCA and the k • p method, which describes the non- 
interacting band structure of pure GaAs. We show that 
nonlocal effects are very important for properly capturing 
the magnitude of Tc, the saturation magnetization, and 
the magnetic anisotropy of this material. In particular, 
we show that the strong cluster anisotropy is responsible 
for the magnetic anisotropy along the cube-edge direction 
and the spin reorientation at low temperature. 
The model Hamiltonian we adopt is 

fl' = i/k.p + Jc5IS(i?/)-J(i?/), (1) 
/ 

where the first term describes the electronic structure 
of the host material (GaAs) in the k • p approximation 
and the second term introduces a magnetic interaction 
between the carrier spin (J) and the local magnetic mo- 
ment (S) of Mn at position i?/. The large magnitude of 
the Mn magnetic moment {S = 5/2) allows us to treat it 
classically. This model is generally accepted to describe 
DMS [3,B i, H liijiik since a mean-field treatment 
of the Hamiltonian^oTllTl [l8| is able to explain many 
physical properties of the system. For the k • p Hamilto- 
nian (-ffk p), we adopt a 4x4 Luttinger-Kohn model de- 
scribing heavy and light hole bands with spin-orbit cou- 
pling, but ignoring the conduction and split-off bands. 
We use the Lut ting er parameters 71 = 6.98, 72 — 2.06, 
and 73 ~ 2.93 [l9|. Biaxial strain is included in Ifk p 
through the strain tensor Exx = £yy = £0 = Aa/a and 
Ezz = 2ci2/cii)eo with the ratio of elastic stiffness 
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FIG. 1: (Color online) (a) Magnetization of Gai_j;Mni,As cal- 
culated by DMFT {Nc = 1) and DCA (iV^ = 16, Nc = 22, 
and Nc = 24) with Mn doping x = 0.05 and hole concen- 
tration p — 0.025. No strain effect is considered, (b) Fer- 
romagnetic transition temperature (Tc) as a function of hole 
concentration (p) with DMFT {Nc = 11 and DCA (iVc = 16) 
for X = 0.05. Experimental results [2a428l | are also shown. 
The Mn concentration in experiments is a;=0.085 (Ku et al), 
a;=0.017-0.09 (Wang et al), and a;=0.05 (S0rensen et al). 



constants C12/C11 = 0.46. Parameter a is the lattice con- 
stant of Gai_i:Mn:rAs, and Aa is the difference between 
the lattice constants of Gai_a;Mna;As and the substrate. 
We use the hydrostatic deformation potential a„ = 1.16 
eV and the shear deformation potential h = — 2.0 eV [l9| . 

In addition to the parameters of the k ■ p Hamiltonian, 
we must determine the value of the exchange coupling 
Jc. It can be obtained from photoemission [20j, infrared 
21, 22 1, and resonant tunneling spectroscopy and 
magneto-transport experiments 24], which give Jc = 0.6- 
1.5 eV. We adopt Jc = 1 eV throughout this Letter. 

Figure[IJa) shows the magnetization per Mn ion of the 
Gai_2;Mn2.As system as a function of temperature with 
DMFT (iVc = 1) and DCA (A^c = 16, A^c = 22, and 
A"c = 24). We chose three fee clusters that are perfect 
according to Belts et al. [2^. The difference between 
DCA and DMFT stems from nonlocal effects, not cap- 
tured in DMFT. The Tc with DCA is far lower than that 
obtained with DMFT, approaching the regime of experi- 
mental values [26I- 28 1 [see Fig.[IJb)]. Another important 
point is the reduction of the saturation magnetization 
at low temperature, consistent with experiments 28-30]. 
This behavior is a product of the rotational frustration 
[J], to be discussed in detail later. This effect also re- 



CD 



90° 
60° 
30° 
0° 
30° 
15° 

0° 
30° 
15° 

0° 



■ ■ '-M-fl-ri-*-^ 

N=\6 - 

N=12 • A 
A.=24 * ^ En = +0.2% ^-f- 



■ Eq = -0.2% 




N=16 m 
N=22 • 
iV^=24 ^ 
1 — 



□an 



(C) 



Welpe?fl/. □ 

Sawicki et al. O 

Masmanidis et al. ^ 

Wang et al. * 



0.2 0.4 0.6 0.: 

TIT^ 



FIG. 2: (Color online) (a) Polar angle {9) of the magneti- 
zation as a function of the normalized temperature {T /Tc) 
for two strain values. Compressive (eo = —0.2%) and tensile 
(eo = -1-0.2%) strain induce in-plane and perpendicular mag- 
netic anisotropy, respectively, (b) Azimuthal angle (</>) of the 
magnetization with respect to the [110} direction with com- 
pressive strain. Experimental results [311 - i34i] are provided in 
(c) to compare with (b). 



duces Tc. The dependence of Tc on hole concentration 
(p) is shown in Fig. [IJb). Tc attains a maximum value 
when hole concentration is half of Mn concentration, con- 
sistent with previous DMFT studies [l^ . 

Next, we studied the magnetic anisotropy of 
Gai_a;Mn2;As. The magnetic anisotropy of this system 
depends on strain, hole concentration, and temperature 
in a complicated manner, but generally it has in-plane 
anisotropy with compressive strain and perpendicular-to- 
plane anisotropy with tensile strain [l|. With compres- 
sive strain, the magnetization changes direction within 
plane from [110] or [110] at high temperature to [100] or 
[010] at low temperature [USi [see Fig. HJc)]. As is 
shown in Fig.[2l^a) and^Jb), the DCA reproduces experi- 
mental results on the dependence of magnetic anisotropy 
on strain and temperature remarkably well. 

While the strain dependence of the magnetic 
anisotropy was explained within the mean-field the- 
ory 171, llSj . the spin reorientation within the plane 
has not been explained yet. In the absence of strain, 
Gai_a;Mn2;As has diagonal magnetic anisotropy within 



the mean-field theory because the heavy holes, which 
dominate at low carrier density, have larger density along 
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FIG. 3; (Color online) (a) DCA results at low temperature 
show that the dimer magnetization (Mdimcr) is prefered to be 
perpendicular to the vector connecting the two Mn ions (R) 
for nearest-neighbor (NN) Mn-dimers. (b) For larger dimers, 
Mdimcr aligns with the total magnetization (Mtot) irrespec- 
tive of R. The dotted and solid arrows represent Mtot and 
Mdimcr, respectively, (c) and (d) Angle between Mdimcr and 
R and between Mdimcr and Mtot, respectively, at r=23.2 K, 
1=0.05, p=0.025, and eo = -0.2%. The left panel is for the 
12 NN dimers [e.g., when the two Mn ions are at (0,0,0) 
and (a/2, a/2, 0)], the middle panel is for the 6 next-nearest- 
neighbor (NNN) dimers [e.g., when the two Mn ions are at 
(0,0,0) and (a, 0,0)1, and the right panel is for the third- 
nearest-neighbor (3' NN) dimers [e.g., when the two Mn ions 
are at (0,0,0) and (a, a/2, a/2)]. 



the diagonal direction in k-space [l3| • With compressive 
strain, it has the anisotropy in [110] or [110]. This ex- 
plains the magnetic anisotropy at high temperature. 

The [100] or [010] anisotropy at low temperature is due 
to the cluster anisotropy originating from the anisotropic 
interaction between neighboring Mn ions. Because Mn 
ions are distributed randomly throughout the system, the 
number of Mn ions within a cluster varies between zero 
and Nc- We call a cluster a monomer (dimer) when it 
includes one Mn ion (two Mn ions). All possible dis- 
tributions of Mn ions are considered effectively in this 
calculation, but at low doping, the magnetic properties 
are dominated by Mn-monomers and Mn-dimers. Since 
there is no cluster anisotropy in monomers, we investi- 
gate magnetization of Mn-dimers in detail. Because of 
translational symmetry, we need to consider only Nc — 1 
dimers. The two Mn ions are nearest-neighbors in 12 
dimers, next-nearest-neighbors in 6 dimers, and third- 
nearest-neighbors in Nc — 19 dimers. Figure [3] shows 
the magnetization direction of each dimer obtained by 
the Monte-Carlo method at low temperature. The mag- 
netization of the nearest-neighbor Mn-dimer is always 
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FIG. 4: (Color online) (a) Magnetic configurations of the 12 
nearest-neighbor (NN) Mn-dimers (Mdimcr) that maximize 
the total magnetization (Mtot) when (j) = Q° (i.e. Mtot is 
along [110] or [110]). The dotted and solid arrows represent 
Mtot and M^mcr, respectively. The numbers below each dia- 
gram indicate the degeneracy of the Mn-dimer configurations. 
The dimer configuration with Mdimcr perpendicular to the 
vector connecting the two Mn ions (R) is energetically fa- 
vored, (b) Same as (a) but when (j> = 45° (i.e. Mtot is along 
[100] or [010]). (c) Maximum value of the average Mdimcr ^s. 



perpendicular to the vector connecting the two Mn ions. 
This cluster anisotropy prevents the magnetic moment 
of some dimers from aligning parallel to the total mag- 
netization and leads to the rotational frustration ^, Iq] . 
When the two Mn ions are farther apart, this anisotropy 
is very weak and magnetization of the dimer aligns par- 
allel to the total magnetization. 

The cluster anisotropy can also explain the enhance- 
ment of Tc up to 260 K in the quasi-two-dimensional 
(5-doped systems 



35]. 



When the Mn ions are within 
one plane, all the nearest-neighbor Mn-dimers can point 
in the same direction (perpendicular-to-plane direction) 
without rotational frustration. This may induce a larger 
saturation magnetization and higher Tc in these systems. 

Notably, due to this cluster anisotropy, the maximum 
total magnetization depends on the magnetization di- 
rection, and this dependence introduces another type of 
magnetic anisotropy. When we assume the magnetic mo- 
ment is perpendicular to the vector connecting the two 
Mn ions and in-plane magnetization, the maximum value 
of the average magnetization of the 12 nearest-neighbor 
dimers is calculated to be 



{Ml 



iincr/ 



COS (j) + sin (j) H- ■\/3 + sin(2(/)) 
hV3-sin(20)l (2) 



for the fee lattice. This becomes maximal when the mag- 
netization is along [100] or [010], as shown in Fig. SJc). 
Since larger magnetization leads to lower magnetic en- 
ergy, it introduces magnetic anisotropy along [100] or 
[010]. We note that this effect becomes unimportant 
at high temperature, where the total magnetization is 
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small. Thus, the spin reorientation from [110] or [110] 
at high temperature to [100] or [010] at low tempera- 
ture is captured within our calculation. This behavior 
arises from the multi-impurity scattering and cannot be 
obtained within the mean-field theory or the DMFT. 

In summary, we investigated the magnetic properties 
of the prototypical DMS system Gai-ajMn^jAs by DCA 
together with the k • p method. We showed that nonlo- 
cal effects, not included in the mean- field theory or the 
DMFT but included in the DCA for Nc> I, are very im- 
portant to quantitatively explain the Tc, the saturation 
magnetization, and the magnetic anisotropy of this ma- 
terial. We find that the edge-direction anisotropy at low 
temperature is due to the cluster anisotropy, and this al- 
lows us to reproduce the spin reorientation in this system 
remarkably well. 
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